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ABSTRACT 

A theoretical model is developed for the interaction of intense femtosecond laser pulses with solid targets on 
the basis of the two-temperature equation of state for an irradiated substance. It allows the description of the 
dynamics of the plasma formation and expansion. Comparison of available experimental data on the amplitude 
and phase of the complex reflection coefficient of aluminum with the simulation results provides new information 
on the transport coefficients and absorption capacity of the strongly coupled Al plasma over a wide range of 
temperatures and pressures. 
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1. INTRODUCTION 

The optical properties of a plasma formed on the surface of an aluminum target subjected to ultrashort laser 
pulses are investigated in interferometric experiments 1-4 at irradiation intensities up to I ~ 10 16 W/cm 2 . In 
the present work, we analyzed data 4 recently obtained with femtosecond pulses of I < 10 W/cm 2 . Such laser 
pulses produce a thin layer of the strongly coupled plasma with an electron temperature up to T ~ 10 eV on 
the initially cold metal target. We developed semiempirical models of optical, transport, and thermodynamic 
properties of Al over the whole range of temperatures realized in a numerical simulation of the laser-metal 
interaction processes. In contrast to previous simulations, 1-3,5 particular attention is focused on the initial stage 
(t ^ 1 ps) of the heating and expansion of the plasma under the conditions of the undeveloped hydrodynamic 
motion of ions. 

Comparison of the numerical simulation results with experimental data on the complex reflection coefficient 4 
makes it possible to determine the free parameters in the formulated expressions for the effective electron collision 
frequency, in contrast to the models discussed elsewhere. 5-7 In contrast to the approach proposed in Ref. 8, this 
study is free of the assumptions that the metal is not melted and the electron temperature is much lower than 
the Fermi energy. This allows us to use the developed models in a broad region of states of the strongly coupled 
plasma formed on the target surface. 

The detailed description of the experimental technique for simultaneous measurements of both the amplitude 
and phase of the complex reflection coefficient can be found in Ref. 4. The source of radiation in the experiments 
of interest 4 was a Cr:forsterite laser system generating femtosecond pulses at a wavelength of Ai = 1240 nm. 9 
The full width at half maximum (FWHM) of the pulses that was measured using an autocorrelator of the 
non-collinear second harmonic was equal to T£ ~ 110 fs in the experiments for the sech 2 envelope shape. The 
time profile of the pulse was measured in a wide range of power using the third-harmonic correlator. The ratio 
(contrast) of the intensity at the pulse maximum to the intensity 1 and 2 ps before the maximum was no less 
than 10 4 and 10 6 respectively. 10 
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For the discussed below experimental data, 4 the target was heated by a p-polarizcd laser pulse at the main 
laser wavelength Ai for the angle of incidence 45°. The spatial distribution of the pumping radiation intensity 
on the target corresponded to the Gaussian with a focus beam diameter of about 70 /zm at a level of exp(— 2). 
A weak s-polarized probe pulse with varying time delay (the second harmonic A2 = 620 nm) was incident 
perpendicularly to the sample surface. 



2. THEORETICAL MODEL 

The self-consistent theoretical model includes the system of electrodynamic equations for describing the ab- 
sorption and reflection of laser radiation, ionization kinetic equations, and one-fluid hydrodynamic equations 
including electron- ion relaxation and electron heat conduction, 11, 12 as well as a new wide-range two-temperature 
equation of state of the irradiated substance. 

For the characteristic time intervals of processes under consideration (t ^ 1 ps), all the typical sizes of 
inhomogencity in the z direction perpendicular to the target surface do not exceed 1 /im. These are much 
smaller than the inhomogeneity sizes along the target surface, which are determined by the size of the focusing 
spot and are equal to tens of microns. Therefore, to analyze the experimental data, 4 we use the one-dimensional 
version of the model developed that takes into account variations of all the quantities only in the z direction and 
the single velocity component V of the quasi-neutral motion of a plasma directed perpendicular to the target 
surface. 

The hydrodynamics equations of continuity for the volume concentration of heavy particles (atoms and ions) 
n a and changes in the substance momentum for the velocity V as well as in the energies of electrons e e and 
heavy particles e l are written in the standard form including thermal ionization in the model of the average ion 
charge Z, 11 

d t n a + d z (n a V)=0, (1) 
[d t + Vd z )V = -p- 1 d z (P e + P i ) + F p , (2) 
n a [d t + vd z y = -P t d z V + Q e \ (3) 

Zn a [d t + Vd z ]e e = -d z q T -Qz + Qib - e e 6 - P e d z V - Q e \ (4) 

[d t + vd z ]z = e. (5) 

Here, P e and Pj are the pressures of electrons and heavy particles respectively, determined by the equation 
of state of the substance; G and Qz are the total rate of the thermal ionization and power density spent on 
ionization respectively, which are calculated using the model 11,12 of the average ion charge, the Lotz formula for 
impact ionization, and the detailed balance principle; p = m,in a is the substance density, m; is the heavy particle 
mass, Z = n e /n a , n e is the concentration of electrons (quasi-neutrality is assumed), F p is the ponderomotive 
force, 

F^-K^n^-dAVzW (6) 

e e|I?|(m e u;i) _1 is the quiver velocity of electrons under the laser irradiation with the slowly varying 
in time amplitude of an electric field E and the frequency uj\, e and m e are the electron charge and mass, 
£ w = 3v / 7r^off(4w) _1 , v e fi is the effective collision frequency of electrons, function 
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arises due to integration of the velocity-dependent cross-section of electron-ion collisions over the Maxwell dis- 
tribution function; 11 Q el — ^ el Zn a {T e —T{) is the electron- ion relaxation energy density, where T e and Ti are 
the temperatures of electrons and heavy particles respectively, and the coefficient j el for lattice temperatures 
Ti < Tmeit (Trndt is the melting temperature) is a constant (j el — 4.93 • 10 10 s _1 for aluminum) and, for higher 
temperatures, is determined by the plasma formula 7" = 3(m e /mi)v e s, where m e is the electron mass and v e g 



is the effective electron collision frequency; qr = K'T e d z T e is the electron thermal flux, where the coefficient 
K 1 for lattice temperatures Tj < T m eit is a constant [K 1 = 4.35 • 10 36 [erg cm s] _1 for aluminum) and, for 
higher temperatures, is determined by the plasma formula, 13 K' = — 128KzZn a /(37rm e i , e s), the factor kz — 0.7 
represents the effect of electron-electron collisions on heat conduction; Qjb is the power density of the inverse 
bremsstrahlung absorption of the energy of the pump laser pulse, 

Qib = (87r)-Wm{ £ (a;i)}|£| 2 , (8) 

where E is the electric-field amplitude of the laser pulse, e is the permittivity of matter, and oj\ is the heating 
(pump) radiation frequency. 

The slowly varying in time amplitude of the electric field strength E = E x e x + E z e z (e x and e z are the 
unit vectors, the plane xOz is determined by the wave vector of incident radiation and the normal to the target 
surface, i.e. direction Oz) for p-polarizcd pump laser pulse was expressed through the amplitude of the magnetic 
field B = Be y as E x = — i/(e(iOi)ki)d z B, E z — —Bsin(9)/s(ui), k\ = uji/c, c is the velocity of light. The wave 
equation (with proper boundary conditions) for the laser magnetic field envelope _B, 14 

d 2 B + fc?[e(wi) - sin 2 (9)}B - d z [ln e(uJi)}d z B = 0, . . 

9,B\ z=Zo =ie{u 1 )k 1 coB(e)[2E L -B]\ s=ao1 BUoo^O, w 



was solved numerically. Here El = \/&kIl/c, II = II if) is the intensity of incident laser pulse, 9 is the angle of 
incidence (the angle between the direction Oz and the wave vector of incident radiation, 9 — 45° for the discussed 
experiments), the point z — zq ^ is at the plasma-vacuum boundary position, and the laser pulse propagates 
from the left (before the laser heats up the target zo = 0, the target is at z 0, i.e. at the right half-space). 

Note that Eq. ((9J describes also a collisionless resonance absorption, but in the discussed range of the pump 
laser intensities, 1^2- 10 14 W/cm 2 , for a high contrast laser pulse of 100 fs duration the characteristic scale 
length of the plasma density does not exceed a few nanometers and the electron plasma temperature is about 
10 eV, and so the laser absorption takes place at over-critical densities in the normal skin heating regime. 15 

The complex reflection coefficient of a weak s-polarized probe laser pulse with frequency L02 at normal inci- 
dence was determined in the linear approximation. For this goal the wave equation 14 with boundary conditions 
for single component of electric field of the probe laser pulse E = E(z)e y , 

d 2 z E + k 2 2 [e(LU 2 ) - sin 2 (9)}E = 0, k 2 = uo 2 /c, , . 

d z E\ z=ZQ = ik 2 cos{6)[2E L -E}\ z=Z0 , E z ^ x = 0, 1 UJ 

was solved numerically with 9 = 0° on the spatially non-uniform profiles of permittivity e(u>2) for different time 
delays. For each time delay the permittivity is determined as described below with the help of numerical solution 
of Eqs. QJ-jSI) that describe the hydrodynamics of a target irradiated by the pump laser pulse determined self- 
consistently by Eq. ([9|). 

It should be noted, that in Eqs. (HJ) and (fl"0|) for subpicosecond laser pulses the transverse (to the normal of 
the target, Oz) gradients of fields and plasma permittivity, which are determined by characteristic scale lengths 
of the order of ten micrometers, can always be omitted in comparison with longitudinal gradients, which are 
determined by the skin-layer depth of a typical size of order of ten nanometers. 

The effective collision frequency z^s over the entire temperature range is determined as the minimum of the 
three values, 

Z4ff = mm{iw,fc'pi,fma X } ) ( n ) 
where v ma t is the effective collision frequency in the metal plasma for T e < Tp = (37r 2 Zn a ) 2 / 3 h 2 /(2m e ), v v \ is 
the collision frequency for the weakly coupled plasma, 16 



v vX = (A/3)V2nZ 2 7i a e i A/ \J m e T e 3 , (12) 

A is the Coulomb logarithm, and i/ maK is the maximum collision frequency determined by the condition 1 that 
the collision mean free path of electrons A pe ~ v e /v c g (v e is the mean electron velocity) is no less than the mean 
distance between ions tq ~ n a -1 / 3 , 

"max = kiWp e , (13) 



where uj pe = \J AnZn a e 2 /m e is the electron plasma frequency and the numerical coefficient ki < 1 is chosen by 
fitting the calculations to the experimental data. 4 Despite the simplicity of such an approach to the determination 
of v e g, it ensures satisfactory accuracy as compared to much more complicated methods. 18 

The effective collision frequency in the metal plasma v me t is determined as 

iW = C e . ph Ti/h + k 2 T 2 /hT F , (14) 
C e . ph = C 00 + C [l - mm{T t /T mclu I} 1 ' 2 ]. (15) 

The first and second terms in Eq. (|14j) represent the contributions from the electron-phonon 19 and electron- 
electron 20 collisions respectively. The constant k 2 is chosen by comparing with the experimental data (also see 
Ref. 8). The first and second terms in Eq. (|15p represent the contributions from the intraband and interband 
transitions 19, 21 respectively. The constant Coo is determined from the data on the static conductivity of the 
metals 21 (for aluminum, Coo — 3.28). The constant Co is chosen so as to ensure the tabulated reflection coefficient 
|r| 2 of the metal under consideration at room temperature. For aluminum, \r\ 2 ~ 0.96 and 0.91, whereas Co — 23 
and 10 for the heating and probe pulses wavelengths Ai = 1.24 pm and A2 = 0.62 fim respectively. When the 
lattice temperature Tj exceeds the melting temperature T mo i t , the band structure of the metal is destroyed 
and the contribution from the interband transitions to the electron-phonon collision frequency vanishes. 22 This 
circumstance is taken into account by a phenomenological dependence on Tj /T mo i t in Eq. I|15p . 

The permittivity of substance is determined by the Drude formula for the metallic plasma 21 for T e ^ T\ = 
0.75Tp and by the formula for a weakly coupled non- degenerate plasma 11 ' 12 ' 16 for T e > T 2 = 1.5T F . In the 
interval T% < T e < T2, a linear interpolation between the Drude formula and plasma formula is used with the 
above-indicated u p \ value. The optical electron mass in the Drude formula is taken as m op t = 1.5m e . 

The thermodynamic characteristics of the condensed phase of the target substance for both the thermal 
equilibrium between ions and electrons and non-equilibrium heating (when T e > Tj) are determined by using a 
new semiempirical equation of state in a wide region of densities and temperatures. In this equation of state, 
the free energy F{p, Tj, T e , Z) is represented as the sum of two terms, F — Fi{p, Tj) + F e {p, T e , Z), determining 
the contributions of heavy particles and electrons respectively. 

The first term, Ft = F c (p) + F a (p,Ti), includes the energy of the interaction between heavy particles and 
electrons at Tj = T e = (F c ) and the contribution from the thermal motion of heavy particles (F a ). The 
dependence of cold energy F c {p) is determined by the procedure described in Ref. 23, which ensures the equality 
of the total pressure in the system to atmospheric pressure at normal density (for aluminum, p = 2.71 g/cm 3 ) 
and room temperature, as well as the agreement with the available data of the shock-wave experiments and 
Thomas-Fermi calculations with quantum and exchange corrections for the high energy densities. The thermal 
contribution of heavy particles to the free energy is given by the expression 24 

3T, , (ei , T,y* 



where a — p/po- To determine the dependence of the characteristic temperature O a = Q a (p), we use the 
interpolation formula 25 



9 (p) = (7 2 / 3 exp 



B 2 + D 2 f B a \na 



(7o a — 2/3) — arctan 



B a \B 2 + D a (\na + D a ) 



(17) 



where 70a is the Griineisen coefficient under normal conditions. The constants T a , B ai and D a are determined 
from the requirement of the optimal description of the experimental data on the thermal expansion and shock 
compressibility of porous samples of the substance. 

The free energy of the electron gas in the metal is given by the expression 

F e (p,T e ,Z) = In (1+^). (18) 

2m,i \ 6Jf / 



Equation (TT8"|) for low and high temperatures is an equation for an ideal degenerate Fermi gas and an ideal 
Boltzmann gas of free electrons respectively. 26 Expressions similar to Eq. (fl"8|) were previously used in the 
equations of state in Refs. 25 and 27. The internal energies of one heavy particle and electron, as well as the 
total pressures for heavy particles and electrons, are expressed in terms of F as e 1 = mi\Fi — Ti(dFi/dTi) p ] and 
e e = mi Z- x [F e - T e (dF e /dT e ) p<z ], as well as P l = p 2 (dF l /dp) Tl and P e = p 2 (dF e j 'dp) Tc ,z respectively. 

3. SIMULATION RESULTS 

Figure [T] shows the results of modeling as well as the experimental data 4 on the absolute value (ri n d) and the 
phase (^ind) of the complex reflection coefficient of aluminum as functions of the maximum intensity I\ of a 
heating laser pulse for various time delays (At = 130, 530, and 930 fs) of the pump and probe pulses. The 
theoretical lines in Fig. [1] are obtained by simulating the heating of the target by the pump laser pulse, and 
also by the calculation of the amplitude and phase of the reflected probe pulse at the doubled frequency using 
described above theoretical model. 

Comparison of the experimental data with the calculations of the absolute value r in( j and phase 'find of the 
complex reflection coefficient of the probe pulse (see Fig. [T]) makes it possible to determine (in the framework of 
the present model) the important properties of the strongly coupled solid-density plasma such as the maximum 
effective frequency of electron momentum relaxation given by Eq. (|13[) and the contribution to the effective 
electron-electron collision frequency specified by Eq. (p~4|) . The best agreement with the experimental data for 
the aluminum targets is reached at ki ~ 0.3, which is close to the theoretical estimate, 17 and at k2 ~ 0.85 
in Eqs. (|13[) and (TTJJ for the effective collision frequencies. In this case, the uncertainty in the choice of the 
coefficients does not exceed 15% including the experimental errors. 

For the indicated parameters, the proposed model well reproduces r; n d and the phase 'find as functions of the 
laser pumping pulse intensity I\ for all time delays At between the pump and probe pulses in the measurements 
with h > 5 • 10 13 W/cm 2 . For fluxes h < 5 • 10 13 W/cm 2 , the model is in agreement with the experimental 
data for At < 500 fs, whereas significant discrepancies exist for r; n d at At > 500 fs. These discrepancies can be 
caused by the formation of a two-component mixture consisting of the low-density plasma and condensed-phase 
fragments on the target surface. 28 A detailed description of this effect is beyond the scope of the model and is 
the subject of further investigations. 

Figure [2] shows the effective collision frequency and the field strength of the laser probe pulse, as well as 
the non-ideality parameter T e i — Ze 2 / (n a ~ 1 ^ 3 T e ) and the degeneration parameter n e A| = (8tt /3)(Tp /T e ) 3 ^ 2 as 
functions of the target depth for the fixed intensity I\ = 6.7 • 10 13 W/cm 2 and the time delay At — 530 fs. Note 
that the average ion charge for the parameters under consideration does not change in the calculation time and 
remains equal to its initial value Z = 3. 

As is seen in Fig. [21 the plasma parameters are strongly inhomogeneous in the skin-layer region, where the 
reflected-signal field of the laser probe pulse is formed. For this reason, despite short duration of the processes 
under consideration, the Fresnel formulas describing the reflection from a homogeneous medium with a stepwise 
boundary are inapplicable to the calculation of the reflection coefficient and its phase. 

Analysis of the simulation results shown in Fig. [5] indicates that the plasma formed on the target surface 
for the experimental parameters under consideration is in a strongly coupled state (T e i > 1) and is strongly 
degenerate (n e A| > 1) in the entire region except for the plasma corona. 

Figure [3] illustrates the formation of a shock wave propagated into the target, which is irradiated by the pump 
laser pulse of I\ = 6.7 • 10 13 W/cm 2 at different time delays. For the times At < 1 ps, both the pressure and 
temperature of electrons are higher than those of heavy particles are. The heavy particle velocity at At = 930 fs 
in the shock wave is of V ~ 0.3 km/s. In the direction opposite to the shock-wave propagation, the plasma is 
expanded with the particle velocity of V ~ 10 km/s. Under the rarefaction, a region of negative pressures takes 
place, see Fig. (3Jb). Such behavior of pressure causes a fragmentation of irradiated target at later stages of the 
considering process. 28 



4. CONCLUSION 



Thus we proposed semiempirical models of optical, transport, and thermodynamic properties of Al over a wide 
range of temperatures and pressures. Those allowed us to have carried out the two-temperature hydrodynamic 
simulation of the processes in the metal under the action of intense femtosecond laser pulses. We emphasize 
that the data on the amplitude and the phase of the reflected field of the probe laser pulse from the experiment 4 
make it possible to acquire important information on the transport properties of the strongly coupled plasma 
formed under the action of the pump laser radiation on the target surface in the subpicosecond time intervals. 
It is shown that the plasma inhomogeneity decisively affects the reflective properties of the target even in such 
short times. 
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Figure 2. Target-depth dependencies calculated for the aluminum plasma characteristics at At — 530 fs and Ji = 
6.7- 10 13 W/cm 2 : (a) normalized (solid line) electric field strength of the probe pulse and (dashed line) effective frequency 
of electron collisions, (b) the parameters of (solid line) non-ideality and (dashed line) degeneration. 



